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We investigate a dynamical mass generation mechanism for the ofT-diagonal gluons and ghosts in 
SU{N) Yang-Mills theories, quantized in the maximal Abelian gauge. Such a mass can be seen as 
«^ ' evidence for the Abelian dominance in that gauge. It originates from the condensation of a mixed 

^ , gluon-ghost operator of mass dimension two, which lowers the vacuum energy. We construct an 

*vj ■ effective potential for this operator by a combined use of the local composite operators technique 

rf^ ' with the algebraic renormalization and we discuss the gauge parameter independence of the results. 

We also show that it is possible to connect the vacuum energy, due to the mass dimension two 
^^ . condensate discussed here, with the non-trivial vacuum energy originating from the condensate 

f^ ' (^m)' which has attracted much attention in the Landau gauge. 

^: 

O ■ PACS numbers: 11.10.Gh,12.38.Lg 



I. INTRODUCTION. 



A widely accepted mechanism to explain color confinement in SU{N) Yang-Mills theories is based on the dual 
Si^ ■ superconductivity picture [l|,l3]i according to which the low energy regime of QCD should be described by an effective 
5h Abelian theory in the presence of magnetic monopoles. These monopoles should condense, giving rise to a string 

formation a la Abrikosov-Nielsen-Olesen. As a result, chromoelectric charges are confined. This mechanism has 
received many confirmations from the lattice community in the so-called Abelian gauges, which are useful in order to 
isolate the effective relevant degrees of freedom at low energy. 

According to the concept of Abelian dominance, the low energy region of QCD can be expressed solely in terms of 
Abelian degrees of freedom p| ■ Lattice confirmations of the Abelian dominance can be found in [J, Q • A particularly 
interesting Abelian gauge is the maximal Abelian gauge (MAG) , introduced in Hi S 13 • Roughly speaking, the MAG 
is obtained by minimizing the square of the norm of the fields corresponding to off-diagonal gluons, i.e. the gluons 
associated with the N{N — 1) off-diagonal generators of SU{N). Doing so, there is a residual U{1)^~^ Abelian gauge 
freedom corresponding to the Cartan subgroup of SU{N). The renormalizability in the continuum of this gauge was 
proven in |3, ISl , at the cost of introducing a quartic ghost interaction. 
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To our knowledge, there is no analytic proof of the Abelian dominance. Nevertheless, an argument that can be 
interpreted as evidence of it, is the fact that the off-diagonal gluons would attain a dynamical mass. At energies below 
the scale set by this mass, the off-diagonal gluons should decouple, and in this way one should end up with an Abelian 
theory at low energies. A lattice study of such an off-diagonal gluon mass reported a value of approximately 1.2GeV 
|lOij. More recently, the off-diagonal gluon propagator was investigated numerically in |ll|. reporting a similar result. 

There have been several efforts to give an analytic description of the mechanism responsible for the dynamical 
generation of the off-diagonal gluon mass. In [i3jll3i ^ certain ghost condensate was used to construct an effective, 
off-diagonal mass. However, in ,14] it was shown that the obtained mass was a tachyonic one, a fact confirmed later 
in [13. Another condensation, namely that of the mixed gluon-ghost operator (^A^A'^" + ac°'c'^) [sj, that could be 
responsible for the off-diagonal mass, was proposed in [I3I. That this operator should condense can be expected on 
the basis of a close analogy existing between the MAG and the renormalizablc nonlinear Curd- Ferrari gauge [l7l Il8l| . 
In fact, it turns out that the mixed gluon-ghost operator can be introduced also in the Curci-Ferrari gauge. A detailed 
analysis of its condensation and of the ensuing dynamical mass generation can be found in |l9l |20| . 

The aim of this paper is to investigate explicitly if the mass dimension two operator (^^4° A'^° -|- ac°'c"') condenses, 
so that a dynamical off-diagonal mass is generated in the MAG. The pathway we intend to follow is based on 
previous research in this direction in other gauges. In [2l|, the local composite operator (LCO) technique was used 
to construct a renormalizable effective potential for the operator A^Af^"^ in the Landau gauge. As a consequence of 
(Af^yl^'^) 7^ 0, the gauge bosons acquired a mass l21j. The fact that gluons in the Landau gauge become massive has 
received confirmations from lattice simulations, see for example |i22] . Recently, th e dy namical mass generation in the 
Landau gauge has been investigated within the Schwinger-Dyson formalism in 1231 l24l| . The condensate (Ai^A^'^) has 
attracted attention from theoretical i23, !2q as well as from the lattice side [23, [23, 123 ■ It was shown by means of 
the algebraic renormalization technique |3^] that the LCO formalism for the condensate (Ai^A'^'^) is renormalizable 
to all orders of perturbation theory jsj . The same formalism was successfully employed to study the condensation of 
{^A^A^^"^ + ac^c^) in the Curci-Ferrari gauge [i^|23. We would hke to note that the Landau gauge corresponds to 
a = 0. Later on, the condensation of A^A^^ was confirmed in the linear covariant gauges J52, J33|, which also possess 
the Landau gauge as a special case. It was proven formally that the vacuum energy does not depend on the gauge 
parameter. However, in practice, a problem occurred due to the mixing of different orders of perturbation theory, 
when solving the gap equation for the condensate. Nevertheless, we have been able to present a way to overcome this 
problem [33|. As a result, it turns out that the non-trivial vacuum energy due to the condensate (^A^A^^^ in the 
Landau gauge coincides with the non-trivial vacuum energy due to the appropriate mass dimension two condensate 
in the hnear covariant gauges, i^A^A^^^), and the Curci-Ferrari gauge, (^\A^A^^^ + ac^c^), since these two classes 
of gauges both have the Landau gauge, a = 0, as a limiting case. 

We would also like to underline that the concept of a gluon mass has already been widely used in a more phe- 
nomenological context since long ago, see e.g. [sj. More recently, a gluon mass of the order of a few hundred MeV 
has been proven to be very useful in describing the radiative decay of heavy quarkonia systems |33| as well as to derive 
estimates of the glueball spectrum |2a| • 

To make this paper self-contained, we will explain all necessary steps in the case of the MAG, and refer to the previ- 
ous papers for more details where appropriate. In section II, we introduce the MAG and discuss its renormalizability 
when the operator (iyl"yl.^° + acf-c"") is introduced in the theory. We briefly review how the effective potential is 
constructed by means of the LCO technique. In section III, we discuss the independence of the vacuum energy from 
the gauge parameter of the MAG. We face the problem of the mixing of different orders in perturbation theory, and 
we provide a solution of it. In section IV, we construct a generalized renormalizable gauge that interpolates between 
the MAG and the Landau gauge. Moreover, we will also show that there exists a generalized renormalizable mass 
dimension two operator that interpolates between the mass dimension two operators of the MAG and of the Landau 
gauge. This can be used to prove that the vacuum energy obtained in the MAG is the same as that of the Landau 
gauge. In section V, we present explicit results, obtained in the case of SU(2) and to the one loop approximation. 
We end with conclusions in section VI. 

II. SU{N) YANG-MILLS THEORIES IN THE MAG. 

Let A^ be the Lie algebra valued connection for the gauge group SU{N), whose generators T^, satisfying 
l^jT^^j^Bj _ jABCj'C ^ g^j-g chosen to be antihermitean and to obey the orthonormality condition Tr (T^^T-^) = 
-TpS^^, with A,B,C =1,..., {N^ - l) . In the case of SU{N), one has Tp = ^. We decompose the gauge field into 
its off-diagonal and diagonal parts, namely 

A^ = A^T"^ ^ A^T" + A'^T', (1) 



where the indices i, j, ... label the TV — 1 generators of the Cartan subalgebra. The remaining N{N — 1) off-diagonal 
generators will be labelled by the indices a, b, .... For further use, we recall the Jacobi identity 

j-ABC rCDE I rADC rCEB , rAEC rCBD _ Q f2) 

from which it can be deduced that 

rabi fbjc I j:abj rbci r\ 

^abc ^bdi i rabd rbic , rabi rbcd r\ /o\ 



The field strength decomposes as 



^M^ = F;^.T^ - F^^T^ + F'^,T\ (4) 



with the off-diagonal and diagonal parts given respectively by 

f;, = DfAt - DfAl + gr'^A^Al , (5) 

f;, = d^Al-d^A'^+gr'^A'I^Al, 

where the covariant derivative D"^^ is defined with respect to the diagonal components Al 



For the Yang-Mills action one obtains 

1 



Df^d^S^'-gr'^A^^ . (6) 



Sym = -iJ d'x {F-^F'^^- + F;,F''-) . (7) 

The so called MAG gauge condition amounts to fixing the value of the covariant derivative, D°^''A'^'', of the off-diagonal 
components by requiring that the functional 

n[A] = (VT)-^ f d'^x {A-I^Af"") , (8) 

attains a minimum with respect to the local gauge transformations. This corresponds to imposing 

DfA"'' = . (9) 

However, this condition being non-linear implies a quartic ghost self-interaction term is required for renormalizability 
purposes. The corresponding gauge fixing term turns out to be |3, 1^1 

Smag =s Jd^x (r (i?;i^A^^ + ^6'') - ^gf ^'^rc'c' - ^gf '^^'c'^Z^r) , (10) 

where a is the MAG gauge parameter and s denotes the nilpotent BRST operator, acting as 

sAl = - {Dfc^ + gf '^"^Ay + gf ^'^A^) , sA^ = - {d^d + gf ^-^A-^) , 

sc" = gf^'^'c^d + ^f^^c^c", sd = ^f'^^c'^c^, 

sc^ = b" , sc' = b' , 

sb" ^ 0, sb' = 0. (11) 

Here c°, c* are the off-diagonal and the diagonal components of the Faddeev- Popov ghost field, while c°, 6" are the off- 
diagonal antighost and Lagrange multiplier. We also observe that the BRST transformations Hll() have been obtained 
by their standard form upon projection on the off-diagonal and diagonal components of the fields. We remark that 
the MAG p()|l can be written in the form 

-'4- (-A^A^"" - -, 
2 ^ 2 



5mag = s7s / d^T -AM'- - ^c^-c'^ , (12) 



with s being the nilpotent anti-BRST transformation, acting as 

sA; = - {Dfc^ + gf '^^'^Allf + gf '^^^A^c^) , sA^ = - (9^c' + gf ^^'A^c") , 

_a ^ gf^bt-b-t j^ Ijabc-b-c^ — < ^ g j iab-a-b ^ 

sc" = -b" + gf^^c^c" + 5/°^V¥ + gf'^'-d'c' , sc' = -b' + gf^^c^lf , 

■sb" = -grH^r ~ gPH^-ff + gP^'tV W ^ -gf^%''c'' . (13) 

It can be checked that s and s anticommute. 
Expression pU|) is easily worked out and yields 

S^MAG = / d'^x (b" (ofM'^ + |6") + c''D''J'D'''"'c'' + gc'T'" (dJ^M^'^) c' + gc''Df l^f-=d^,.c^d^^ 

- agfH'^-d'd - g^r'^f'^^'c^c'^AiA'"' ~ -gf^H'^c^c" - ^g^f^b^Jcdi^a^b 



_, .__ ^ , , c c c c 

_ ^g2j^abcjadi-b^^d^t _ ag2jabcjade-b-c^dA ^ ^^^^ 
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We note that a = does in fact correspond to the "real" MAG condition, given by eq.©- However, one cannot set 
a = from the beginning since this would lead to a nonrenormalizable gauge. Some of the terms proportional to a 
would reappear due to radiative corrections, even if a = 0. See, for example, |37l |. For our purposes, this means that 
we have to keep a general throughout and leave to the end the analysis of the limit a ^ 0, to recover condition ^. 
The MAG condition allows for a residual local [/(1)^~^ invariance with respect to the diagonal subgroup. In order 
to have a complete quantization of the theory, one has to fix this Abelian gauge freedom by means of a suitable further 
gauge condition on the diagonal c omp onents A^ of the gauge field. A common choice for the Abelian gauge fixing, 
also adopted in the lattice papers |10l [ll| , is the Landau gauge, given by 

5diag = s I d^x c'd^A^" = / d^x {b'df.A^^' + c'd^ (9^c' + gf ''''A^^c'')) , (15) 

where c*, 6* are the diagonal antighost and Lagrange multiplier. 

A. Ward identities for the MAG. 

In order to write down a suitable set of Ward identities, we first introduce external fields f7'^% fl^"", U, L° coupled 
to the BRST nonlinear variations of the fields, namely 



^o.-xt — 



f d^x {-f}^" {D'^J'c'' + gt'^'Alc" + gr'"Alc') - fl^"' {d^,c' + gP^^Alc}") 



with 



L" [gn'd'c' + ^/'^""cV") + L'2. p'"'d'c!'\ , (16) 

sn"" ^ sfl'"' = , (17) 

sL"" = sV = . 

Moreover, in order to discuss the renormalizability of the gluon-ghost operator 

1 

2" 
we introduce it in the starting action by means of a BRST doublet of external sources (J, A) 

sA = J, sJ = 0, (19) 

so that 

^AI.A^- + ac-cA + C^ 



Omag = t;^^''" + "^"C^ ' (18) 



^LCO = s I d'^x [X[ -A'^^A^- + ac'^c'^ ) + C— ) (20) 



r2 



= I d^x {j{ ^A'^^A^'' + ac^cA + C Y - aA&'^c" 



^^t^a (^jjab^b _^ gJab^J^b^^^^ ^ ^^-a fj aU ^b ^^ _^ gj 



bi ^b ^i I y f abc b c 







^a,^ 


■^a.r 


&"•' 


A 


J 




j^a,r 


dimension 
ghost number 


1 






1 


2 

-1 


2 



2 
-1 


2 



3 
-1 


4 

-2 



TABLE I: Dimension and ghost number. 



where C, is the LCO parameter accounting for the divergences present in the vacuum correlator (OMAG(a^)C'MAG (?/))> 
which are proportional to J^. Therefore, the complete action 



S = ■S'ym + S'mag + Sd 



iag ^ *^cxt 



'5'lco 



(21) 



is BRST invariant 



sE = . (22) 

As noticed in |lfill38l|. the gluon-ghost mass operator defined in ea. l|18|l is BRST invariant on-shell. 

In Table I, the dimension and ghost number of all the fields and sources are listed. We are now ready to write down 
the Ward identities needed to discuss the renormalizability of the model. It turns out that the complete action S is 
constrained by 

• the Slavnov- Taylor identity 



with 



sm 



d^x 



5(S) = , 
ST. 5T 5T 5T ST. ST 5T ST 



(5f7^" (5 A" Snt^"- SAl SL°- (5c" (5L* S& 



fe" 



ST 

s¥ 



,,ST 
Sc' 



s\ 



• the diagonal ghost equation 



where 



G'T 
S 



cl > 



gr'c' 



ST 



and 



A! 



-d^c' + gT'^n'^'A'' - d^n"' - gr'"L''d 



Notice that expression (|27|l . being linear in the quantum fields, is a classical breaking. 
• the diagonal gauge fixing and anti-ghost equations 



— - 8 At'' 



• the integrated A-equation 



Sc"- snt"^ 



ST 



= 



ST 



d^xi -rr + c"—- 



0, 



SX Sb'^ 

expressing in a functional form the on-shell BRST invariance of the gluon-ghost operator Omag- 
• the diagonal U{1)^^^ Ward identity 

W'T ^ -d%^ , 

with 

ST 






L"- 



SLt 



This identity follows from the diagonal ghost equation H25I) and the Slavnov- Taylor identity H23|l . 
In order to find the foregoing Ward identities, use has been made of the Jacobi identities Q. 



(23) 
(24) 

(25) 
(26) 
(27) 

(28) 



(29) 

(30) 
(31) 



B. Algebraic characterization of the most general local counterterm. 

We mention that all the classical Ward identities of the previous section can be extended to all orders of perturbation 
theory without encountering anomalies. In principle, this can be proven by means of the algebraic setup of |30| | and of 
the general results on the BRST cohomology of gauge theories [33 ■ It can be understood in a simple way by observing 
that pure Yang-Mills theory in the MAG can be regularized in a gauge invariant way by employing dimensional 
regularization. 

In order to characterize the most general invariant counterterm which can be freely added to all orders of perturba- 
tion theory, we perturb the classical action S by adding an arbitrary integrated local polynomial E'^"""' in the fields 
and external sources of dimension bounded by four and with zero ghost number, and we require that the perturbed 
action (E + jyE™"'^*) satisfies the same Ward identities as E to the first order in the perturbation parameter 77, i.e., 

5(E + 77E'=°""') = +0(77^) , 
e'(E + 77E^™"') = A^i +0(r]2). 



n,, \ /n , nroiintx « ^/ 2\ 



y'rf4a;(^^-fc''^)(E + r;E— ') = + 0{ii') , 

W(E-^77E™""') = -d^b' +0{Tf). (32) 

This amounts to imposing the following conditions on E'^"""* 

where By, denotes the nilpotent linearized operator 

SsSs = , (34) 



6v = / d 



j4. 



X 



f 6Y. 5 5Y. 5 (5E 5 5Y. 5 6Y. 5 



and 



5Y. 5 SJ: S SY. S 5 ^ S S 



/-»zy^count r\ 

"^^ = 



(35) 



5U 
^From the conditions (|33|l and (|36|l . it turns out that the most general invariant counterterm can be written as 

j^count ^ z^ /■ ^4^ (f;^f''''" + f;^f'^'") + SsA-^ , (37) 

where A~^ is an integrated local polynomial with ghost number —1, given by 

+ aiX (^A'^^A'"' + acfcA + ^XA'^^A'"'' + iaa^Xlfc'' + ^^j") . (38) 



We see thus that S'^"""* contains six free independent parameters, namely qq, ai, 03, 05, a^ and 013. These parameters 
can be reabsorbed by means of a muhiphcative renormahzation of the gauge coupling constant g, of the gauge and 
LCO parameters a, C, and of the fields = (A^°, A''', c'', c", c\ c\ b"", }f) and sources $ = (f7^", ^2^^ L"-, L\ X, J), 
according to 



with 



and 



with 



S(g, a, C, 0, $) + ^S™""* = S(ffo, ^o, Co, ^0, <J>o) + 0{r,') , (39) 

3o = Zgg , ao == ^^a , Co == ^cC > 

A^" = Zy^^'^° , A^^ = Zg-M''* , (40) 

eg = Zy^c-^ , eg = Zl^^r , (41) 

c^ = Zl/'e , cl = Z-V2c« ^ (42) 

6S = ZgZy^^y^^'' , 6^ = Zg6' , (43) 

0^° ^ z^^/^z-iz-1/2^'^'' , n^' = z-1/2 f^M^ ^ (44) 

Lg = Z;1Z-i/2z-i/2l« , LI = Zg-^Z-iL' , (45) 

Jq = Z^a Z^ J = Z~ZcJ , 

Ao = Z^:z-i/2a = ZgZy^A , (46) 



Zg = 1-7] 



Oo 

2 ' 



/2a5 

Z„ = 1 + ?y h ao - 2a6 

V a 

Z(^ = 1 + rj {ai3 + 2ao — 2ai — 2aQ) 

1/2 _ -, , ^fa^ 

2 



^i^ = l + ^(^+«3) , 



zy2 ^ i + i(„g_ai), 

Zy2 = i + |(ag + ai). (47) 

In particular, from ea. (|46|l . one sees that the renormahzation of the source J, and thus of the composite operator 
CmaGj can be expressed in terms of the renormahzation of gauge coupling constant and of the diagonal ghost. This 
property follows from the diagonal ghost equation (|25|l and from the the i nteg rated A-equation (|29|l . In particular, 
for the anomalous dimension of the gluon-ghost operator Omag, we obtain }19J 

70MAo(ff') = ^^j^ log {Z'gZ,) = 2 (^^ - 7c. (5^)) , (48) 

with 

1A9') = A^l^l^^c/'- (49) 



C. The effective potential. 

We present here the main steps in the construction of the effective potential for a local composite operator. A more 
detailed account of the LCO formalism can be found in jifl, llj . 

To obtain the effective potential for the condensate (Omag)! we set the sources J7' J7" i", U and A to zero and 
consider the renormalized generating functional 

exp(-iW(J)) ^ f[D^]cxpiS{J) , 

S{J) = Sym + Smag + 5diag + S-count + / d^x (z J J (^^ZaA^A^'' + Zo^ZcOc'^A + (C + K)^\ (50) 

where ^p denotes the relevant fields and S'count is the usual counterterm contribution, i.e. the part without the 
composite operator. The quantity 5C, is the counterterm accounting for the divergences proportional to J^. Using 
dimensional regularization throughout with the convention that d = 4 — e, one has the following identification 

CoJo' = M"^(C + '5C)J'. (51) 

The functional Vy(J) obeys the renormahzation group equation (RGE) 



O r\ r\ r K f) \ 



(52) 



where 



^From eq.jnU, one finds 



7a(g^) = ^i—\■aa^^l—\nZ^'^, 

vig^O = A^I^C- (53) 



77(5^ C) = 270MAO (g'K + Sig', a) , (54) 



%^a)= (e + 27oMAo(5')-/3(5')^-«7a(5')^)<SC. (55) 



with 



Up to now, the LCO parameter C is still an arbitrary coupling. As explained in |40ll41| . simply setting C = would 
give rise to an inhomogeneous RGE for W{J) 



(56) 

and a non-linear RGE for the associated effective action T for the composite operator Omag- Furthermore, mul- 
tiplicative renormalizability is lost and by varying the value of 5C,^ minima of the effective action can change into 
maxima or can get lost. However, C, can be made such a function of g^ and a so that, if g^ runs according to P{g'^) 
and a according to 7a(5^), C{g^,ct:) will run according to its RGE l|^ . This is accomphshed by setting ( equal to 
the solution of the differential equation 

/?(5')g^ +«7a(5',a)— j ag^,a) = 2^o^,^i9')ag' ,a) + 5{g\a) . (57) 

Doing so, yV{J) obeys the homogeneous renormahzation group equation 

'm|^ + Pig')^ + "^"(5')^ - 70mao(3') / d'xJ^^ W{J) = . (58) 

To lighten the notation, we will drop the renormalization factors from now on. One will notice that there are terms 
quadratic in the source J present in yV{J), obscuring the usual energy interpretation. This can be cured by removing 



the terms proportional to J^ in the action to get a generating functional that is linear in the source, a goal easily 
achieved by inserting the following unity, 



1 = jj I [Da] exp iJd'^xU^f^^ Omag ~ C,J 



2> 



(59) 



with N the appropriate normalization factor, in ea. H5U|) to arrive at the Lagrangian 



L[A,,,o) - --Fl.F^-'^ - -Jl.F^''' + /:mag + U^.^ - Tu ^ -^9<^Omag - -, (^mag)' , (60) 



-F^^F*^""- - -^FJi^^F'^"' + Cmag + -Cdiag - ^-5^ + -j^gcrOuAG - ^ 
while 



exp(-iW(J)) = f[Dip]expiSa{J) , (61) 

Sa{J) = I d^x(c{A^,cj) + J-\ . (62) 

^From eas. ()5(J|l and H(jlfl . one has the following simple relation 

5W{J) 



5J 



= - (Omag) = -(-), (63) 

j=o \ 9 ! 



meaning that the condensate (Omag) is directly related to the expectation value of the field a^ evaluated with the 
action Sa = J d'^x£{Afj_, a). As it is obvious from eq-JHOJ), (a) 7^ is sufficient to have a tree level dynamical mass for 
the off-diagonal fields. At lowest order (i.e. tree level), one finds 



off — diag. 
^'^gluon 



-gLf ^- - J-'-^- (64) 



Meanwhile, the diagonal degrees of freedom remain massless. This could have been established already from the local 
[7(1)^-1 Ward identity ^. 

III. GAUGE PARAMETER INDEPENDENCE OF THE VACUUM ENERGY. 

We begin this section with a few remarks on the determination of Cifj^i ct). From explicit calculations in perturbation 
theory, it will become clear fSql that the RGE functions showing up in the differential equation H57I) look like 

70mag(5^) = lQ{a)g'^ + 7i{a)g* + ■ ■ ■ , 

laig^) = ao{a)g'^ + ai{a)g'^ -\ , 

%^a) = So{a) + Siia)g^ + ■ ■ ■ . (65) 

As such, ea. (|57|) can be solved by expanding C(.9^i ct) in a- Laurent series in g^, 

^^g^^a)^^ + <:i{a)+C2{a)9' + --- • (66) 

g 

More precisely, for the first coefficients (0, Ci of the expression l|t)t)|l . one obtains 

a/- 

2/3oCo + aao— — = 270C0 + Sq , 
oa 



2f3iCo + aao-^ \- aai— = 270C1 + 271C0 + ^i 



(67) 
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Notice that, in order to construct the n-loop effective potential, knowledge of the (n + l)-loop RGE functions is 
needed. 

The effective potential calculated with the Lagrangian 1)60(1 will explicitly depend on the gauge parameter a. The 
question arises concerning the vacuum energy -Bvac, (i-G- the effective potential evaluated at its minimum); will it be 
independent of the choice of a? Also, as it can be seen from the equations 1(671) . each Ci(a) is determined through a 
first order differential equation in a. Firstly, one has to solve for Co(a)- This will introduce one arbitrary integration 
constant Cq. Using the obtained value for Co(Q!), one can consequently solve the first order differential equation for 
Ci(a). This will introduce a second integration constant Ci, etc. In principle, it is possible that these arbitrary 
constants influence the vacuum energy, which would represent an unpleasant feature. Notice that the differential 
equations in a for the Ci are due to the running of a in eg. ((57(1 . encoded in the renormalization group function 70(9^). 
Assume that we would have already shown that -Evac does not depend on the choice of a. If we then set a = a* , with 
a* a fixed point of the RGE for a at the considered order of perturbation theory, then equation 1(57(1 determining ^ 
simplifies to 



f3i9')^C{9', a) = 2jo^,,{g')ag', a*) + 5{g\ a*) 



(68) 



since 



7a (3^ 



)a 



0. 



(69) 



This will lead to simple algebraic equations for the Q{a*). Hence, no integration constants will enter the final result 
for the vacuum energy for a = a* , and since E^ac does not depend on a, E'vac will never depend on the integration 
constants, even when calculated for a general a. Hence, we can put them equal to zero from the beginning for 
simplicity. 

Summarizing, two questions remain. Firstly, we should prove that the value of a will not influence the obtained 
value for -Evac- Secondly, we should show that there exists a flxed point a*. We postpone the discussion concerning 
the second question to the next section, giving a positive answer to the first one. In order to do so, let us reconsider 
the generating functional 1(61(1. We have the following identification, ignoring the overall normalization factors 



exp(-iW(J)) = f[Dtp]cxpiSaiJ) = ^ f[DtpDa]expi S{J) + f 



2CU 



d^x\ -—{ --Omag-CJ 



(70) 



where S{J) and S'cr(J) are given respectively by eg. lfSUI) . and ea. (|H^ . Obviously, 

2 



i^/[^"]^"P 



^ I d^xi-^f ^ ^Omag-CJ 



d 
da 



(71) 



so that 



dWjJ) 
da 



s I d xs \ —Cc" 



terms ex J 



(72) 



j=o 



which follows directly from 



dS{J) _ /■ 4 /I ^ „ 

; — ss a X \ -c c 

da J \2 



terms ex J 



(73) 



We see that the first term in the right hand side of ((73(1 is an exact BRST variation. As such, its vacuum expectation 
value vanishes. This is the usual argument to prove the gauge parameter independence in the BRST framework 
[3(1 |. Note that no local operator O, with sO = Cmag, exists. Furthermore, extending the action of the BRST 
transformation on the a-field by 



fabc—abc 



sa = gsOuAG = -A''''D'l^''c'' + ab^c" - ag f'""c''-S' d - -gr'"'-&'c''c'' 



(74) 



one can easily check that 



f I d'^xC{A^,(j) =0 



(75) 
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so that we have a BRST invariant cr-action. Thus, when we consider the vacuum, corresponding to J = 0, only the 
BRST exact term in ea. (|72|l survives. The effective action F is related to WiJ) through a Legendre transformation 



- ) = -W(J) - 



d'yJ{y) 



o-(2/) 



The effective potential V{(j) is then defined as 



Via) d^x^Tl- 



(76) 



(77) 



From 



dV{a) 
da 



(78) 



(f) 



it follows that 



and hence, we derive from eas. H76(l and (|77|l that 



^"W 



F = -J 



J = 



d'^x 



i»'<^> 



j=o 



Thus, due to ea. (|7^ . 



da 



Via) 



(79) 



(80) 



(81) 



(82) 



We conclude that the vacuum energy E'vac should be independent from the gauge parameter a. 

A completely analogous derivation was obtained in the case of the linear gauge (231 • Nevertheless, in spite of the 
previous argument, explicit results in that case showed that i^vac did depend on a. In |23| it was argued that this 
apparent disagreement was due to a mixing of different orders of perturbation theory. Let us explain this with a 
simple example. Let us first notice that a key argument in the previous analysis is that the source J = vanishes at 
the end of the calculations. In practice, J = is achieved by solving the gap equation (f7H|l . Moreover, in a power 
series expansion in the coupling constant, the derivative of the effective potential with respect to a will look like 



{vo + v,g^ + 0(/)) a , 
where we assume that we work up to order g^. The corresponding gap equation reads 



Due to eas. ((77)l and ifT?^ . one also has 
Imposing the gap equation l|84() leads to 



Vo + vig^ + Oig^) = 



J = 5 («o + «iff' + 0(/)) . 



J^g{0 + Oig^))a 



(83) 



(84) 



(85) 



However, as it can be immediately checked from expression (|7U|) . there are several terms proportional to J in the 
right-hand side of eq.jT^Jl. For instance, one of them is given by ^J^- Since 



5C ^5Col 
da da g"^ 



da 



Oig' 



(87) 



12 



we find 



d^j2 
da 



^Co 2 , /^^Co„ , f^Cl 2 \ 2 , ^/ 4n \ 2 



Squaring tlie gap equation H84|l. 



leads to 



_^„^-^ _|2.„„.H-3^„„^|.' + 0(.-)). 



i;^ + 2vivog'^ + 0(/) = 






da 



da 



(88) 



(89) 



(90) 



We see that, if one consistently works to the first order, terms such as ^ J^ do not equal zero, although J = to that 
order. Terms like those on the right-hand side of ea. (|90|l are cancelled by terms which are formally of higher order, 
requiring thus a mixing of different orders of perturbation theory. Of course, this problem would not have occurred 
if we were be able to compute the effective potential up to infinite order, an evidently hopeless task. Nevertheless, 
in [331 we succeeded in finding a suitable modification of the LCO formalism in order to circumvent this problem 
and obtaining a well defined gauge independent vacuum energy i^vac, without the need of working at infinite order. 
Instead of the action (|50|l . let us consider the following action 



S{J) — Sym + •S'mag + Sdi 



diag 



d^X 



JTig^,a)OMAG + ^J'H9^a)J^ 



where, for the moment, T{(j^ ^ a) is an arbitrary function of a of the form 

T{g\ a) = 1 + /o(a)g2 + /^(a)/ + o(/) 



(91) 



(92) 



and J is now the source. The generating functional becomes 



exp(-i>V(J)) = / [1)0] expiS'(J) . 



(93) 



Taking the functional derivative of Y<}{^J) with respect to J, we obtain 

(5W(J) 



<5J 



= -T{g\a){p^KG) . 



.7=0 



Once more, we insert unity via 



1 - - / {Do\ exp 



to arrive at the following Lagrangian 



i \ d^x\ -^ 



1 / a 



2C V5-^(5^«) 



OMAG-C^^(5',a) 



L{A^,,a) = --F^^F^''" - ^i^^^i^^''' + /:mag + ^diag - ^ 2jf2(„2 



(94) 



(95) 



2g-T\g\a)C, + ^W, «)C^''''"^« " ^ (0MAG)^(96) 



^From the generating functional 



exp(-i>V(J)) = [D<j>]expiSa{J) 



S^{J) = / d^x /:(A^,5) + J- 



it follows that 



SWjJ) 
SJ 



j=o 



{a)^gT{g\a){0 



MAG/ 



(97) 
(98) 

(99) 
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The renormalizability of the action (|62(l imphes that the action (|98|l wiU be renormahzable too. Notice indeed that 
both actions are connected through the transformation 

The tree level off-diagonal masses are now provided by 

off — diag. 
"^gluon 

K"'- = ^/«T^' (101) 

while the vacuum configuration is determined by solving the gap equation 

dV(a 




da 



, (102) 



with Via) the effective potential. Minimizing Vja) will lead to a vacuum energy £'vac(«) which will depend on a and 
the hitherto undetermined functions fi{a) |56|. We will determine those functions fi{a) by requiring that i5vac(Q^) is 
a-independent. More precisely, one has 

dEys_c 

— ; — - = =^ first order differential equations in a for /i(a) . (103) 

da 

Of course, in order to be able to determine the fi{a), we need an initial value for the vacuum energy -Evac- This 
corresponds to initial conditions for the fi{a). In the case of the linear gauges, to fix the initial condition we employed 
the Landau gauge [S^], a choice which would also be possible in case of the Curci- Ferrari gauges, since the Landau 
gauge belongs to these classes of gauges. This choice of the Landau gauge can be motivated by observing that 
the integrated operator / d'^xA^A^^^ has a gauge invariant meaning in the Landau gauge, due to the transversality 
condition d^^A^^^ = 0, namely 

[VT]-^ min I d^x {A'^)^ {A^'^)^ = f d^xiA^Af"^) in the Landau gauge , (104) 

with the operator on the left hand side of ea. (|104|l being gauge invariant. Moreover, the Landau gauge is also an 
all-order fixed point of the RGE for the gauge parameter in case of the linear and Curci- Ferrari gauges. At first 
glance, it could seem that it is not possible anymore to make use of the Landau gauge as initial condition in the case 
of the MAG, since the Landau gauge does not belong to the class of gauges we are currently considering. Fortunately, 
we shall be able to prove that we can use the Landau gauge as initial condition for the MAG too. This will be the 
content of the next section. 

Before turning our attention to this task, it is worth noticing that, if one would work up to infinite order, the 
expressions (|OT|l and (|^ can be transformed exactly into those of (jSnil . respectively (|^ by means of ea. HlUU|l and 
its associated transformation 

a = T{g^,a)a, (105) 

so that the effective potentials V{a) and V{a) are exactly the same at infinite order, and as such will give rise to the 
same, gauge parameter independent, vacuum energy. 

IV. INTERPOLATING BETWEEN THE MAG AND THE LANDAU GAUGE. 

In this section we shall introduce a generalized renormahzable gauge which interpolates between the MAG and the 
Landau gauge. This will provide a connection between these two gauges, allowing us to use the Landau gauge as 
initial condition. An exaniple of such a generalized gauge, interpolating between the Landau and the Coulomb gauge 
was already presented in |42j |. Moreover, we must realize that in the present case, we must also interpolate between 
the composite operator ^^Af^A^^ of the Landau gauge and the gluon-ghost operator Omag of the MAG. Although 
this seems to be a highly complicated assignment, there is an elegant way to treat it. 
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Consider again the SU{N) Yang-Mills action with the MAG gauge fixing (|12() . For the residual Abelian gauge 
freedom, we impose 

- Kgf''^A'^^Ai'''{h'' -gP^''c''c^)^Kgf''''A^\D'l^c'')t ~^ng^r^'r'''^A'^^A^'''c'^-^) , (106) 

where k is an additional gauge parameter. The gauge fixing H1Q6|I can be rewritten as a BRST exact expression 



•-"diag — d X 

Next, we will introduce the following generalized mass dimension two operator. 



(107) 



l^a^^a + |. 



O = -A'^^A'"' + -A^A*"' + a-ec" , (108) 



by means of 



^Lco = s^d^x (^XO + C^^ (109) 

- KXc'df.A^'' + ngf^^XAI^A^^'c^ 

with (J, A) a BRST doublet of external sources, 

sX = J, sJ = 0. (110) 

As in the case of the gluon-ghost operator H18|l , the generalized operator of eq. H108() turns out to be BRST invariant 
on-shell, a property which can again be expressed in a functional way, see ea. (|121f) . 
Let us take a closer look at the action 

S = Sym + S'mAG + 'S'diag + 'S'lCO + '5'ext • (HI) 

The external source part of the action, 5extj is the same as given in ea. p6|l . 

Also, it can be noticed that, for k — > 0, the generalized local composite operator O of ea. (|108|) reduces to the 
composite operator Omag of the MAG, while the diagonal gauge fixing (|107|l reduces to the Abelian Landau gauge 
(|15|l . Said otherwise, for k ^ 0, the action E' of ea. (|lll|l reduces to the one we are actually interested in and which 
we have discussed in the previous sections. 

Another special case is k ^ 1, a — > 0. Then the gauge fixing terms of E' are 



5mag + ^^iag = / d^xs i-A^d^c^) = J d^x {c^d^^Dfc^ + b^d^A 



(112) 



which is nothing else than the Landau gauge. At the same time, we also have 

1 



lim O ^ -A^A^'^ , (113) 



which is the pure gluon mass operator of the Landau gauge |2lL l31( . 

^From 31], we already know that the Landau gauge with the inclusion of the operator A^A^^^ is renormalizable 
to all orders of perturbation theory. On the other hand, in section II, we have proven the renormalizability for k = 0. 
Before we continue our argument, let us first prove the renormalizability of E' for general a and k 7^ 0. The complete 
action E', as given in eq.l^Hj, is BRST invariant 

sE' = , (114) 

and obeys the following identities 
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• The Slavnov- Taylor identity, provided by 

5(S') = / d'^x 



4^ / 5E' (5E' (5E' 5E' ,5E' (5E' (5E' (SE' 
,,<5E' ^,5E' ^<5E'\ ^ 



(115) 



• The integrated diagonal ghost equation 



where 



and 

a classical breaking. 
• The diagonal anti-ghost equation 



and 



0'E' = A 



cl ; 



G' 



(Tx 



5 



gr'c- 



ct^x [gf'''''n^"'Al - gf^^'L^c^ + K\^^A^'' 









(116) 



(117) 



(118) 



(119) 



(120) 



• The integrated generalized A-equation 



d^x 



5 a S r S' 



E' ==0 



(121) 



expressing the on-shell BRST invariance of the operator O of ea. pU8|l . 

Also in this case, these Ward identities extend to the quantum level. Accordingly, the most general local countcrtcrm 
j^/count ]^yg^ obey the following constraints 





Bi 


yi /count 


= 0, 


/J zy /count 


= 0, 


^ v^ /count 

6¥^ 


= 0, 




y>/count 


= 0, 


c 


+ C 

[sx Sb-\ 


y /count 


= 0. 



d^x 

where Bs' denotes the nilpotent, Bs'B^' ~ 0, linearized operator 

(5E' S (5E' S 



By 



d^x 



ST.' S ST.' 6 ST,' S 

+ 



g^at, Sj^a Sj^a g^at, g^,,', gM gM g^^,^ gj^a g^a 



ST/_J_ ST/__5_ ST^_S_ a±_ ,i±_ j± 
5c'' SL'^ SL"- S& S& SU Sc" 5'^ SX 



(122) 



(123) 
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^From general results on BRST cohomology |33, we know that the most general, local counterterm can be written as 

(124) 



where A"-"^ is an integrated local polynomial of ghost number —1 and dimension 4, given by 



a'^niA"" + a'^niAi"' + a'^L^c" + a'^Vc' + a'g {d^c") A^° + a'lg f^'c'' AlA^' 



+ a'^ (5^c') A'" + a'^gr^lftc^ + ^agr'"rc''c' + ^agr''^rc''c^ + aa'^^b^^r 



The constraints (|122|l lead to the relations 

oy = a'li = a'ls = , 



(125) 



06=02, 

^14 = 2a-|^Q + 03 , 
a'i2 = a'i- a's , 

a'l = 4 + K(ai3 -4) , 
/ / 

02 = 04 = and thus 0^3 = . 



(126) 



Summarizing 

A — 1 / j4_ r I oa Alia I I ra a I —a j-^ab Atih I rabi—a ai Aub f rabi—a—b i 

= lax \ai\l^A^ + a^g^L c ~ a^c U^ A'^ — na^gj c A^A^ — an^agj c c c 



——agj c c c + aa^QO c + a^A I — A A^ + ac c I —AAA'^ + ZaiQaAc c -\ — — (AJ 



2 



(127) 



In comparison with the case of the MAG, we see that J]"^"^"* also contains six free independent parameters, namely 
Og, a'l, a3, 05, a'lQ and a'l^, despite the fact that the action E' contains the extra gauge parameter k. These parameters 
can be reabsorbed by a suitable multiplicative renormalization of the gauge coupling constant g, of the gauge and 
LCO parameters a, k, (, and of the fields cj) = {A^^"-, A^^i, c°, c", c*, c*, 6°, V) and sources $ = (fJ'"", fi^*, L°, L*, A, 
J), according to 



where 



E'(5, a, K, C, 0, $) + ?7E'^°"'^* = S'(5o, ao, ^o, Co, 0o, $o) + 0(r,2) , 



50 = Zgg , ao = ^aOf , Co = ^cC , ^0 = ^0 ^^a ^g ^'* > 



(128) 



^7 = z)['^A'"' , vi;;' = z-^Af"' 



(129) 



Cq — Z/c C , 



-a _ ^l/Zr^a 
Cq — Z/j, c , 






6S = Zgzy2zi/2j,a 



r!^'' = z^'/'z-iz-i/'r!'^'^ 



&0 = Zgh' , 



j^M^ ^ ^-1/2 ^^^^ 



(130) 
(131) 
(132) 
(133) 
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Lg - Z-'Z-y-'Z^y-'L'^ , LI = Z-'Z-'V , (134) 

(135) 



with 



7-17-1/ 


'2^-l/2^a ^ 


4 = Z;iZ-iL^ 


J^ = z 


gZcJ , Ao = 


= z,zl/'x , 


z, = 


1-4, 




Zc. = 


1 + ri (24 + aj, - 


-<) , 


Zc = 


1 + V {a'm + 2ao 


+ 24 - 24) , 


zT = 


l + Jf + aC 


). 


z'J' = 


i-|K + 4), 




z'J' = 


i + |K-4)- 





(136) 

We see thus that the additional gauge parameter k does not renormaUze in an independent way. Furthermore, from 
ea. (|135|l . we notice that the relation (|48|l is generalized to the operator O, i.e. 

^0(9') - ~2 [^ - 1A9')) . (137) 

Summarizing, we have constructed a renormalizable gauge that is labelled by a couple of parameters (a, k). It allows 
us to introduce a generalized composite operator O, given by ea. pU8|l . which embodies the local operator A'j^A^^ 
of the Landau gauge as well as the operator Omag of the MAG. To construct the effective potential, one sets all 
sources equal to zero, except J, and introduces unity to remove the J^ terms. A completely analogous argument 
as the one given in section III allows to conclude that the minimum value of V{cr), thus £'vac, will be independent 
of a and k, essentially because the derivative with respect to a as well as with respect to k is BRST exact, up to 
terms in the source J. This independence of a and k is again only assured at infinite order in perturbation theory, so 
we can generalize the construction, proposed in section III, by making the function J^ of ea. H92|l also dependent on 
K. The foregoing analysis is sufficient to make sure that we can use the Landau gauge result for E^ac as the initial 
condition for the vacuum energy of the MAG. Moreover, we are now even in the position to answer the question 
about the existence of a fixed point of the RGE for the gauge parameter a, which was necessary to certify that no 
arbitrary constants would enter the results for E^ac- We already mentioned that the Landau gauge, i.e. the case 
{a, k) = (0, 1), is a renormalizable model |3l|, i.e. the Landau gauge is stable against radiative corrections. This can 
be reexpressed by saying that {a, k) = (0, 1) is a fixed point of the RGE for the gauge parameters, and this to all 
orders of perturbation theory. 

V. NUMERICAL RESULTS FOR SU{2). 

After a quite lengthy formal construction of the LCO formalism in the case of the MAG, we are now ready to 
present explicit results. In this paper, we will restrict ourselves to the evaluation of the one-loop effective potential 
in the case of SU{2). As renormalization scheme, we adopt the modified minimal substraction scheme (MS). Let us 
give here, for further use, the values of the one-loop anomalous dimensions of the relevant fields and couplings in the 
case of SU{2). In our conventions, one has p3 . l44Ll45| 

7c^(5') = (-3 - a) ^ -f 0(5'*) , (138) 



167r2 
3 a) 167r2 



7.(5^) = {-2a+l~-)-f^ + 0{g^) , (139) 



while 



.2 ./22 g^ 



Pig') = -eg' - 2 ( -^ ) + O(g^) , (140) 
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and exploiting the relation (|48|l 



70mag(5^) = ( -T^2a 



26 
3" 



IGtt^ 



+ 0(/) , 



(141) 



a result consistent with that of [4J| . 

The reader will notice that we have given only the 1-loop values of the anomalous dimensions, despite the fact that 
we have announced that one needs (n + l)-loop knowledge of the RGE functions to determine the n-loop potential. 
As we shall see soon, the introduction of the function J^{g^ , a) and the use of the Landau gauge as initial condition 
allow us to determine the 1-loop results we are interested in, from the 1-loop RGE functions only. 

Let us first determine the counterterm S(^. For the generating functional W( J), wc find at 1-loop [53| 



W{J) = I d''x[^{C + 6C)Y) + ^Indet [S"'' (O^ 



aJ) 



- In det 
2 



gab 



(92 + J) g^^ _ ( 1 _ i ) d^d. 



and employing 

In det 
with 



S'^" (52 + J) g^. 



dudy 



5"" [{d - l)trln {d^ + j) + trln {d^ + aj)] , 



5"" = N{N - 1) = 2 for TV = 2 , 
one can calculate the divergent part of ea. (|142|l . 



W(J) 



d^x 



'SC 



J' 



3 1 



167r2 



1 2,2! 
-a J - 



2 72 



167r2 



87r2 



raV 



Consequently, 



^<-i(°' 



3)i + 0(/) 



Next, we can compute the RGE function S{g^, a) from ea. (|55|l . obtaining 

^'2-3 



S{g',a) 



a 



87r2 



Oig' 



Having determined this, we are ready to calculate Co- The differential equation (|67|l is solved by 

(o(a) ^a+{9-4a + 3a^)Co , 



(142) 



(143) 



(144) 



(145) 



(146) 



(147) 



(148) 



with Co an integration constant. As already explained in the previous sections, we can consistently put Co — 0. 
Here, we have written it explicitly to illustrate that, if a would coincide with the 1-loop fixed point of the RGE for 
the gauge parameter, the part proportional to Co in ea. (|148|l would drop. Indeed, the equations 9 — 4a -I- 3a^ ~ 
and —2a + | — ^ = 0, stemming from ea. H139() . are the same. Moreover, we also notice that this equation has only 
complex valued solutions. Therefore, it is even more important to have made the connection between the MAG and 
the Landau gauge by embedding them in a bigger class of gauges, since then we have the fixed point, even at all 
orders. In what follows, it is understood that Co = ct- 

We now have all the ingredients to construct the 1-loop effective potential Vi (?) , 



^Vi-f 

2Co V V 


2/0 + ^)5')+* In det 


\ Co 


)] 


In det 

2 


. VV Co/ V 


1\ V 
1 dndu 


J 



or, after renormalization 

;2 



V.(.).^(l-(2/o.^),^ 



g'^' 



32^2 ^2 



In 






J_gW^ 
32^2 ^2 



In 



gacr 



(149) 



(150) 
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We did not explicitly write the divergences and counterterms in ea. (|151|) . since by construction we know that the 
formalism is renormalizable, so they would have cancelled amongst each other. This can be checked explicitly by using 
the unity of (|5(^ with counterterms included. It can also be checked explicitly that Viia) obeys the renormalization 
group 



^^^^(^) 



+ terms of higher order 



by using the RGE functions (|138|l - (|141|l and the fact that the anomalous dimension of a is given by 



tA/) 



Pig' 



+ lOMAoia^) + 1^ 



dlTiT{g^,a) 



which is immediately verifiable from ea.H99|l. 

We now search for the vacuum configuration by minimizing Vi(ct) with respect to a. We will put Jl^ 
possibly large logarithms, and find two solutions of the gap equation 



(151) 



(152) 



^ to exclude 

Co 



dVi 

da 







"«0 



^ ^ ( 1 - ( 2/o 



Ci 

Co 



1-9 



9 o- 



167r2 Co^ 



g o- 



1 g'^a^a 



32^2 ^2 



IBtt^ 



So 



1 3 

ma 

2 



1 g^a'^a 



327r2 



So 



= 



namely 



g'N 



167r2 



2Co 



N=2 



167r2 (2/oCo + Ci) + «2 Ina - a2 ^ 1 



(153) 

(154) 
(155) 



The quantity y is the relevant expansion parameter, and should be sufficiently small to have a sensible expansion. The 
value for (a) corresponding to ea. l|155l) can be extracted from the 1-loop coupling constant 



g'm 



1 



/?oln^ 



(156) 



The first solution (|154|l corresponds to the usual, perturbative vacuum (i?vac = 0), while ea. H155l) gives rise to a 
dynamically favoured vacuum with energy 



£'vi 



off — diao; 



647r^ 



fo 2\ ( off— diagX 

(3 - a ) (^mgi^^„ s j 



gluon 
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Ea. (|157() is obtained upon substitution of ea. (|153l) into ea. (|150() . From ea. (|157() . we notice that at the 1-loop ap- 
proximation, a^ < 3 must be fulfilled in order to have -Evac < 0. In principle, the unknown function /o(a) can be 
determined by solving the differential equation 



fl-C'vac i\ ^ n I off-diag^ 
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da 
dy dfo 
dfo da 



= 




(159) 



with initial condition £^vac('3) — E, 



Landa 



However, to solve ea. (|159|) knowledge of (^i is needed. Since we are not 



interested in fo{a) itself, but rather in the value of the vacuum energy Svac, the off-diagonal mass m 



off- 
gluon 



and the 



expansion parameter y, there is a more direct way to proceed, without having to solve the eq.H159|) . Let us first give 
the Landau gauge value for i^vac in the case N = 2, which can be easily obtained from |2lLl47|. 
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FIG. 1: The off-diagonal giuon (fat line) and ghost mass (thin line) in function of a. Masses are in units of Aj^. 



Since the construction is such that _Evac(a) = -^'vac : ^^ '^'^'^ equally well solve 






(161) 



which gives the lowest order mass 
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and hence 
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(163) 



The result (|162|l can be used to determine y. From ca. (|158|l one easily finds 

36 



y 



187 + 66 In 5^3^ 



(164) 



We see thus that, for the information we are currently interested in, we do not need explicit knowledge of Ci and /q. 
We want to remark that, if ^i were known, the value for y obtained in ea. H164|l can be used to determine /o from 
ea. (ll55|l . This is a nice feature, since the possibly difficult differential equation H159|) never needs to be solved in this 
fashion. In Fig. 1, we have plotted the off-diagonal gluon mass (|162ll and ghost mass (|163ll for < a < -\/3- We 
notice that the masses grow to oo for increasing a, while the expansion parameter y drops to zero, as it is clear from 
Fig. 2. The relative smallness of y means that our perturbative analysis should give qualitatively meaningful results. 
Before we come to the conclusions, let us consider the limit a — *■ 0, corresponding to the 
One finds 



"real" MAG DfA^"^ = 0. 
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(165) 



VI. DISCUSSION AND CONCLUSION. 



The aim of this paper was to give analytic evidence, as expressed by ea. (|165|) . of the dynamical mass generation for 
off-diagonal gluons in Yang-Mills theory quantized in the maximal Abelian gauge. This mass can be seen as support 
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FIG. 2: The expansion parameter y as a function of a. 



for the Abelian dominance 13,0.13 in that gauge. This result is in quahtative agreement with the lattice version of the 



p^+m^ 



MAG, were such a mass was also reported llO|,lll|. The off-diagonal lattice gluon propagator could be fitted by 
which is in correspondence with the tree level propagator we find. We have been able to prove the existence of the 
off-diagonal mass by investigating the condensation of a mass dimension two operator, namely {^A°^A^°- + ac°c°). It 
was shown how a meaningful, renormalizable effective potential for this local composite operator can be constructed. 
By evaluating this potential explicitly at 1-loop order in the case of SU{2), the formation of the condensate is favoured 
since it lowers the vacuum energy. The latter docs not depend on the choice of the gauge parameter a, at least if 
one would work to infinite order in perturbation theory. We have explained in short the problem at finite order and 
discussed a way to overcome it. Moreover, we have been able to interpolate between the Landau gauge and the MAG 
by unifying them in a larger class of renormalizable gauges. This observation was used to prove that the vacuum 
energy of Yang-Mills theory in the MAG due to its mass dimension two condensate should be the same as the vacuum 
energy of Yang-Mills theory in the Landau gauge with the much explored condensate (Ai^A''^). It is worth noticing 
that all the gauges, where a dimension two condensate provides a dynamical gluon mass parameter, such as the 
Landau gauge J2j], the Curci- Ferrari gauges z(J|, the linear gauges 66\ and the MAG, can be connected to each other, 
either directly (e.g. Landau-MAG) or via the Landau gauge (e.g. MAG and linear gauges). This also implies that, if 
(AfJ^yl^'^) 7^ in the Landau gauge, the analogous condensates in the other gauges cannot vanish either. Then the 
question arises if this correspondence between different gauges could be stretched further to for instance the Coulomb 
gauge, where the possibility of a condensate (^AfA^'^'^ was already advocated some time ago in [48j. However, it is 
worth remarking that this might be a more complicated task, since the Goulomb gauge is not a covariant gauge fixing, 
and as such its analysis within the algebraic remormalization framework |3Cl| is not straightforward. 

Needless to say, the present work is far from being complete. First of all, an explicit calculation at two loop order 
and for general gauge group SU{N) would be interesting. We also limited our computations to the tree level order. In 
principle, one should evaluate the off-diagonal gluon polarization in order to get further information on the structure 
of the propagator. A first step in this direction was taken in the case of the Landau gauge in [ig. It is unknown 
what will happen at higher orders in the MAG, but it is likely that the external momentum Q^ will enter through 
loop corrections and influence the possible position of a pole in the propagator. The ghost condensation, that was 
first investigated in |12 . Ili^ as a possible mechanism behind the off-diagonal mass, and later on was shown to be 
tachyonic [hi. Il5j| . could enter this polarization too. This would require a more complete treatment of the ghost 
condensation in the MAG, along the lines of 49], where these condensates were considered in more detail in the case 
of the Curci- Ferrari and Landau gauge. Another issue which deserves attention is the behaviour of the diagonal gluon. 
In [ll|, it was found that the diagonal gluon propagator also contain a mass parameters, with rn J^^^ ~ i^giuon'**^' 
while in |10| the diagonal gluon was reported to behave like a light or massless particle. For completeness, we remind 
that these lattice simulations were both performed in the case of SU{2). We want to remark that a condensation 
of the composite diagonal operator A* A*** cannot occur within our approach, since this is forbidden by the diagonal 
local t/(l)^~^ Ward identity (|3U|) . In principle, one could add an extra source term like ^pA^ A^^ , but it does not 
seem possible to prove the renormalizability of this operator in the MAG. This might be consistent with the result of 
[lll| , since the diagonal gluon propagator could not be fitted with a Yukawa propagator ■i}„^'i , in contrast with the 
off-diagonal gluon propagator which could be fitted with il^-i ■ This could mean that the diagonal mass parameter 
is of a different nature compared to the off-diagonal one. A possible speculation is that it might have to do with 
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Gribov copies, since a fit i^^i did work for the diagonal propagator |11|. 

Our analysis of the MAG condensate was also restricted to the purely perturbative level. One could imagine 
calculating in a certain non-trivial background. The vacuum energy calculated in one gauge should still be the same 
as the one calculated in the other gauge. In this context, and keeping in mind that monopole condensation is an 
essential ingredient of the dual superconductor picture, it might be worth noticing that the role of (A^ Af^^) as an 
order parameter for monopole condensation was investigated in the Landau gauge by the authors of [2y|, based on a 
similar observation in compact QED 25]. We note that an off-diagonal gluon mass can serve as a starting point to 
derive low energy (dual) Abelian models for Yang-Mills theories, see for example [50, |5ll |53 • 

Let us conclude with a few considerations on the issues of the degrees of freedom and of the unitarity when the gluons 
attain a dynamical mass, as a consequence of a nonvanishing dimension two condensate (Cmag) = {^A°i^A^^°' + ac°c°). 
One possible way to look at the degrees of freedom associated to a given field is through its propagator. From the pole 
of the propagator one gets information about the mass of the field, while from its residue one learns about polarization 
states. However, the propagation of the field has to occur in some vacuum. In other words, the kind of vacuum in 
which the field propagates has to be supplemented. In our case, this task is achieved by the LCO Lagrangian, eq.ljSOJj 
i.e. 



l^a^j,,.a _ 1^.^^,.,: ^ ^^^^ ^ ^^^^^ _ _^ ^ ^^gaOuAG " ^ (O.^.r.^' 



which allows one to take into account the effects related to having a nontrivial vacuum corresponding to the nonvan- 
ishing dimension two condensate (Omag); as expressed by the identity 

(a) ^ g (Omag) ■ (167) 

That this is the preferred vacuum follows from the observation that the vacuum energy is lowered by the condensate 
(Omag)- Expanding thus around {a) 9^ 0, a dynamical tree level mass rn°i^^^^'^^' for the off-diagonal gluons is 
generated in the gauge fixed Lagrangian Hl(j6|) . namely 



off-diag. _ / g(^/ /-|fin^, 

^gluon - y ^— • (168) 

Therefore, in the condensed vacuum, (a) 7^ 0, the Lagrangian Hlt)6|) accounts for off-diagonal massive gluons. However, 
we emphasize that this dynamical mass parameter occurs as the result of a particular condensate. It is not a free 
parameter of the gauge fixed theory, its value being determined by a gap equation. Concerning now the unitarity 
of the resulting theory, it should be remarked that, due to confinement, gluons and quarks are only to be called 
physical at a very high energy scale Q^, where they behave almost freely and asymptotic states can be related to 
them, thanks to asymptotic freedom. At very high energies, our dynamically massive action might be unitary: a 
renormalization group improvement could induce quantum corrections such that the mass parameter runs to zero for 
Q^ ^ 00. Otherwise said, the corrections induced by this dynamical mass on the scattering amplitudes are expected 
to become less and less important as the energy of the process increases, so that the amplitudes of the massless case 
are in fact recovered. Such a scenario would be analogous to the behaviour of the dynamical mass parameter discussed 
by Cornwall in [53 ■ For very high Q^, one does indeed expect that perturbative Yang-Mills theory with massless 
gluons, having two physical degrees of freedom, describes the physical spectrum and that non-pcrturbative corrections 
are absent. 

To decide if our resulting theory is unitary at smaller Q^, one should know how to take into account the effects of 
confinement, which now cannot be neglected. This would amount to knowing how to construct out of our Lagrangian 
(|166|l the low energy spectrum of the theory, which is believed to be given by colorless bound states of gluons and 
quarks as, for instance, mesons, baryons and glueballs. This task is far beyond our capabilities. At intermediate Q^, 
what we can state is that this dynamical mass parametrizes the behaviour of the Greens function of the gluon. As a 
result of quantum effects, i.e. the condensation of the mass dimension two operator, a pole appears in the off-diagonal 
gluon propagator at the tree level. Including higher order effects will alter the propagators behaviour as well as 
the location of the pole at physical Q^ (i.e. Q^ < 0). In the case of the Landau gauge, higher order calculations 
showed that the condensate remains stable, and hence a nonzero mass parameter will remain, see |2lll47j . This mass 
parameter will describe the behaviour of the Greens function at Euclidean Q^. The presence of a mass parameter 
does however not necessarily entail the presence of a pole in the propagator at physical Q^. Using lattice simulations 
of the Euclidean propagator, a mass parameter is found also by fitting at Euclidean Q^ > 0, but no pole and related 
massive particle is implied. Analogously, one should not conclude from our calculations that the gluon is a massive, 
physical particle and that unitarity is violated. 
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